





NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3792 


A THEORETICAL ANALYSIS OF HEAT TRANSFER 
IN REGIONS OF SEPARATED FLOW 
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SUMMARY 


The flow field analyzed consists of a thin, constant pressure, 
viscous mixing layer separated from a solid surface by an' enclosed region 
of low- velocity air ("dead air" ) . The law of conservation of energy is 
employed to relate calculated conditions within the separated mixing 
layer to the rate of heat transfer at the solid surface. This physical 
model is applied to laminar separations in compressible flow for various 
Prandtl numbers, including consideration of the case where air is injected 
into the separated region. Application to turbulent separations is made 
for a Prandtl number of unity in low- speed flow without injection. All 
calculations are for the case of zero boundary-layer thickness at the 
position of separation. 

For laminar separations the differential equations for viscous flow 
at arbitrary Mach number are solved for the enthalpy and velocity profiles 
within the thin layer where mixing with dead air takes place . Results 
are presented in tabular form for Prandtl numbers between 0.1 and 10. 

The rate of heat transfer .to a separated laminar region in air (Pr = 0.72) 
is calculated to be O.56 of that to a corresponding attached laminar 
boundary layer having the same constant pressure. Injection of gas into 
the separated region is calculated to have a powerful effect in reducing 
the rate of heat transfer to the wall. It is calculated that a moderate 
quantity of gas injection reduces to zero the heat transfer in a laminar 
separated flow. 


INTRODUCTION 


Separated laminar layers (sometimes called "free layers") have long 
been known to be quite unstable in subsonic flow. A rule of thumb based 
on experiments of von Doenhoff (ref . l) is that a separated laminar layer 
in subsonic flow will remain laminar downstream of separation only for a 
run of about 50,000 in Reynolds number (based on length along separated 
layer),, after which transition occurs. Because of this low stability 
relative to attached laminar flows (which remain laminar to Reynolds num- 
bers the order of 10 6 ), the special type of separated flow wherein a 
viscous layer remains laminar not only downstream of separation but also 
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for a short distance downstream of reattachment has occurred only rather 
rarely in the past. As a result , the "pure laminar" type of separated 
flow, which is the main subject for analysis in this report, formerly has 
been regarded primarily as a laboratory curiosity rather than a practical 
phenomenon . 


Some interest recently has been stimulated in the pure laminar type 
of separated flow, inasmuch as free laminar layers appear to be surpris- 
ingly stable at hypersonic Mach numbers . It was observed in reference 2 
that the stability of separated laminar layers increases markedly with an 
increase in Mach number until, at Mach numbers near A, they can become 
for certain configurations almost as stable as an attached laminar bound- 
ary layer (laminar to about 1 million Reynolds number) . This strong 
increase in stability makes separated laminar flows of practical interest, 
especially at hypersonic Mach numbers . The actual extent of this practi- 
cality, however, is not known at present inasmuch as the available data 
are meager above Mach numbers of about 4. The effect of cold- wall condi- 
tions on the stability of separated laminar flows is not known, nor is 
the effect of gas injection on stability known. 

In view of the trend of increasing stability with increasing Mach 
numbers, it is pertinent to inquire into the heat-transfer characteristics 
of separated flows . It might be expected intuitively that the rate of 
heat transfer from a region of separated laminar flow would be smaller 
than that from a comparable attached laminar boundary layer. This expec- 
tation, in fact, led Stalder and Nielsen (ref. 3) to conduct experiments 
with separated flow induced by a probe in front of a blunt-nosed object. 
Instead of finding a lower heat-transfer rate when separation was induced, 
they observed a doubling in the heat-transfer rate. In retrospect, how- 
ever, it has been concluded from a detailed examination of their original 
shadowgraphs that the separated flow induced by the probe was of the 
"transitional" type (transition between separation and reattachment), 
whereas the attached flow to which they made comparison was that of a 
completely laminar boundary layer. They compared a partly turbulent 
separated flow with a fully laminar attached flow. Hence, the question 
remains open as to whether a pure laminar separated flow has a higher or 
lower rate of heat transfer than a comparable attached laminar flow. 

The main object of the present paper is to compare by means of theo- 
retical calculation the heat-transf er rate in a separated laminar region 
with that of a corresponding attached laminar boundary layer. Secondary 
objectives are to consider the case of gas injection, and to apply the 
theoretical ideas also to turbulent separation. In a previous paper 
(ref. 4), the velocity profiles within a separated laminar mixing layer 
were calculated under the assumption that the specific heat was constant. 
In the present paper, calculations are made of the enthalpy profiles as 
well as the rate of heat transfer. The theory of reference 4 is extended 
to include the case where the specific heat is a function of temperature, 
since this case corresponds to conditions encountered in hypersonic 
flight . 
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NOTATION 

width of two-dimensional flow 
specific heat at constant pressure 

constant of proportionality between viscosity and temperature 

constant defined by equation (24) 

functions defined by equations (23) 

function defined by equation (22) 

functions defined by equations (26) 

enthalpy per unit mass 

functions defined by equations (48) and (49) 
length of separated mixing layer 
mass flux (density x velocity x area) 

Mach number 
static pressure 

total heat flux to wall between separation and reattachment points 

Q w 

average heat flux to wall per unit area, 

0 surface area 

recovery factor based on enthalpy (see eq« (42)) 
radius of base of cone 

Reynolds number based on conditions at outer edge of mixing layer 
and on length l 

gas constant 

Prandtl number 

-static temperature 

local velocity components in x and y directions, respectively 
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x,y coordinates parallel and perpendicular, respectively, to direction 
of flow along dividing streamline (\jr = 0) within mixing layer 

p mass density 

(i viscosity 

t stream function defined by equation (4) 

(j. 

v kinematic viscosity, — 

A coefficient of heat conductivity 

7 ratio of specific heats 

£ dimensionless mass-flow variable defined by equation (10) 

T] function defined by equation (15) 

functions defined by equation (32) 

Subscripts 

a active degrees of freedom of molecule (translational and rotational) 

d dead-air region 

e outer edge of viscous layer 

int internal degrees of freedom within molecule 

i injected mass flow into separated region 

o dividing stream line „ 

t total conditions 

w wall 

* dimensionless quantities defined by equation (7) 

bZ boundary layer 

Superscript 


ordinary differentiation 
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ANALYSIS 
Physical Model 


The type of separated flow analyzed is depicted schematically in 
figure 1. The pressure is regarded as essentially constant over the 
length l of the two-dimensional laminar mixing layer. Also, in accord- 
ance with the theoretical description advanced in reference 2, the region 
of recompression through the reattachment zone is assumed to he of small 
dimensions compared to l . This physical picture yields results consist- 
ent with the experiments of reference 2, in the sense that the dead-air 
pressure calculated according to such a physical picture agrees well with 
measurements for the case when reattachment occurs on an inclined plate . 
This establishes some confidence in the mechanism postulated, at least as 
a reasonable approximation. 


The differential equations of momentum, continuity, and energy for 
viscous flow within the relatively thin, constant -pres sure, laminar 
mixing layer are as follows: 


<3u c)u 

pu + pv 

ox ay 



(1) 


5(pu) d(pv) _ 
dx 3y 



( 2 ) 

(3) 


The conditions at the outer edge of the mixing layer provide two 
boundary conditions: u(x,oo) = u e and h(x,<») = h e , where u e and h e are 

constants since the pressure is constant. The condition that dead air 
borders the low-velocity portion of the mixing layer provides two more 
boundary conditions: u(x,-oo) = 0 and h(x,-o°) = h w . These latter two 

boundary conditions require some explanation since the mixing layer 
cannot extend to -00 due to the presence of the wall. Actually, the use 
of —oo in the boundary conditions means that the mixing-layer character- 
istics are determined as though the wall were an infinite distance away. 
It is assumed, in effect, then, that the presence of the wall distorts 
the mixing-layer profiles only to a minor degree in the region where the 
velocity is low. This will be a good assumption if the dimensions of 
the dead-air region are large compared to the mixing-layer thickness. 
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The above system of differential equations and boundary conditions 
is to be solved subject to the following additional assumptions: 

(i) Pr = constant 

(ii) p = pRT 

(iii) \i/\i e - CT/T e 

(iv) mixing-layer thickness at separation is 0 (or negligible 
compared to l ) 


Assumption (i) is a common one. Assumption (ii), taken together with 
thermodynamics theory, implies that the enthalpy is a function of tempera- 
ture only. This allows the specific heat Cp to vary with- the tempera- 
ture as would be the case in hypersonic flight. Assumption (iii), 
introduced in reference involves a constant C which is selected so 


p re f CT re f 

that — = . if T-r^f is selected as the wall temperature, and 

M ' e Te /tJ7 t + s 

if Sutherland’s equation for viscosity is employed, then C = /~ - 

. T e Tv + b 


Assumption (iv) is made to simplify the mathematical computations, since 
it implies that the profiles of velocity and enthalpy are similar at every 
station within the thin mixing layer. Brief discussion is presented later 
of how the analysis would be modified to consider the case of finite thick- 
ness of boundary layer at separation. 


Transf ormation of Equations 

The solution of equations (l) to ( 3 ) proceeds exactly the same as in 
references 4- and 5 by slightly modifying the von Mises transformation 
(ref. 6) as generalized to compressible flow by von Karman and Tsien 
(ref. 7). The first step is to employ the stream function t as an 
independent variable, where 


pu = p. 


pv = -Pe - 


since this variable automatically satisfies the continuity equation (2) . 
Transformation is made to (x*,t*) coordinates through the transformation 
formulae 

(kl ■ £ (,) 


/V\ = + 

W/„ p s Wy, W/i 
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The following dimensionless variables are introduced: 
u* = u/ue h* & h/h e 
x* — x/z — ^/•JvqMqIQ 

P* = p/p e P* = = CT/T e = CT* 


(7) 


When substituted into equations (l) and (3), the transformation formulae 
and dimensionless variables yield two basic equations 



It is noted that neither the Mach number nor the ratio of specific heats 
7 appears in these equations . Only the ratio Ue 2 /he appears . In the 
special case where Te is low enough (less than about 700° R for air) 
so that no vibrational energy is excited within molecules, then 
u e 2 /h e = (7 e - l)Mg 2 . In the following, however, the form containing 
ue 2 /h e is retained throughout. 

A relatively simple solution of equations (8) and ( 9 ) is made pos- 
sible since the momentum equation (eq. (8)) is independent of the energy 
equation (eq. (9))* The energy equation, however, is not independent of 
the momentum equation and, hence, the momentum equation must be solved 
first . 


Solution of Momentum Equation 


A solution of the momentum equation already has been obtained in 
reference This solution corresponds to similar velocity profiles at 
each station along the mixing layer. The solution of reference 4 is 
given in terms of a single variable 

i = %/j^ Uo) 

which reduces equation (8) to the ordinary differential equation 
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Solution of this nonlinear differential equation, subject to the boundary- 
conditions u*(oo) = 1 and u*(-oo) = 0, was obtained through the numerical 
procedure described in reference The function u^.(^) is tabulated in 
table I . 

A relation between £ and the geometric coordinates (x,y) is 
required later in order to convert the u*((;) function into physical 
velocity profiles u*(x,y) . The appropriate relationship is derived 
from the stream-function differential 


£H dy _ 111 dx = di|r = 
P e p e 


pv 


& 


a£ +f|t 

dx 


dx 


by considering a fixed value of x 


pu dy = p q^VqUqIC d£ - p e ^ v e u e xC d£ 


and then integrating 

y / u e = i r r T 

2 "/ v e xC 2 J o P*U # J o * 2u* 


(12) 


(13) 


It is seen that even though the momentum equation is completely solved' by 
the function u*(£), actual velocity profiles cannot be calculated until 
the temperature function T^(^) is calculated, as this function together 
with the u*(£) function determines the relation between £ and y. 


The variable £ is clearly a dimensionless mass-flow variable 
(see eq. (12)). Inasmuch as the mass flux m drawn out from the dead 
air by the scavenging action of the mixing layer up to station x must 
be finite, 


/ o r£=° 

pu dy = — = ^PeUgPgXC / d£ = finite 

£=£d 


- - £ d \l Pe u eP-e x ^ 




where is the value of £ at u* = 0. Table I shows that =-l- 2 33» 

I, 5 ins' : i.ve, though not necessary to subsequent analysis, to 
digress empo./:,:. . d in order to illustrate the relation between the 
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mass-flow variable £ and the more familiar Blasius variable q 
where 



(15) 


If a function f (n) is employed as the dependent variable instead of the 
function u* (5), & different differential equation is obtained in place 
of equation (ll) . By setting 

f (n) - 5 ( l6 > 


so that 


i f(n) 

the differential equation becomes 


f ,n + 

which is the familiar Blasius equation for incompressible , viscous, 
lam inar flow over a flat plate. The boundary conditions to which equa- 
tion (l8) is subjected, however, are not the same as the familiar Blasius 
boundary conditions. The boundary conditions for a separated laminar 
mixing layer are 


1 d£ 

2 dq 


= u J 


ff” = 0 


(17) 


(18) 


f’(oo) = 2 f (0) = 0 f f (-oo) = 0 


The third of these boundary conditions differs from the corresponding 
boundary conditions in Blasius’ problem, which would be f’(o) = 0. The 
relation between the function f(q) and the function u*(£) is given by 
equations (l6) , (17)> and the relationship f n (q) = Mu^du^/d^ . 


Solution of Energy Equation 


Inasmuch as the variable £ is pertinent to the solution of the 
momentum equation, the energy equation (9) also is transformed to (^,x^.) 
coordinates . 
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u* 


3% 


+ 


Pr L+ 

2 d l ) 


5h * 

hi 



( 19 ) 


This equation will he solved subject to the boundary condition 
h(-«») = h w = constant , which implies that the temperature of the dead air 
is constant just before it is drawn into the mixing layer. It is not 
known whether a constant temperature of this dead air can be realized in 
practice, inasmuch as appropriate experiments have not yet been conducted. 
The case tu^ - constant, though, is relatively simple since the term 
involving o/bx* can be dropped from equation ( 19 ) • If subsequent experi- 
ments show a variable "wall* 1 enthalpy h v to be essential, then equa- 
tion ( 19 ) with the d/dx* term included would have to be solved in a 
manner analogous to the solution of reference 5* In the present analysis, 
however, it is assumed that dh^/dx* =0, so that 


1 


a 

d £ 2 


(Pr| 


OUjAdh* 



( 20 ) 


This is a linear, second-order, ordinary differential equation, and is 
amenable to general solution. By noting from equation (ll) that 
u*u* M + OV? + £u x . , /2 = 0, the general solution to equation (20) can 
be written as 


MO = 1 + c 


00 

1 f (**VV) 

J £ 


Pr 


iL 

2u v 


Pru# 


kh< 


f {ku 


■* u * 


)Pr G(u lr (21) 


where C x is an arbitrary constant to be evaluated from boundary condi- 
tions, and where 


G(0 = J (4u*u*') 


2 -Pr d£ 
2u J1 


( 22 ) 


Actually, the lower limit could, with equal generality, be an arbitrary 
number; zero is a convenient limit for reasons which will later appear. 
For convenience, the following two functions are introduced which will 
appear often in subsequent equations: 

r°° Pr d£ 

F x (?) s / (4u*u»') — 

\ * > 

F 2 (0 (^*u # ’) Pr G(£) ~ 


( 23 ) 
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The constant Ci can now be written as 


h w 



1 ' 2h e 

Fl(«d) 


(2k) 


where F x (£ d ) and F 2 (£ d ) depend only on the Prandtl number. Their values 
have been calculated from equations ( 23 ) and the known solution for u*(£) 
as follows: 



^■fe- 1 + GS-0 8l(e)+ 5£%( £ > (25) 


where the normalized functions gi(£) and g 2 (£) are defined as 

' 

Si(C) = Fi(£)/Fi(£ d ) 

- 

g a (£) -F a (S) - g x (£)F a (5 d ) 

- 

in order to have the following boundary values 


(26) 


gi(£ d ) = 1 '> giW = 0 

S 2 (^d) = 0 > 8 2 (") = 0 

The functions g x (£) and g 2 (£) have been calculated from equations (23) 
and (26) using the numerical solution for u*(£)* These functions have 
been determined for Prandtl numbers of 0.1, 0.25, 0.5, 0.72, 1.0, 1.5, 

2.0, and 10.0. Results are presented in table II. This table enables 
the enthalpy distribution to be calculated readily as a function of £ . 

The Pr range of interest for gas flows is between about 0.25 and 1.0; 
the other values would correspond to liquids (when the density is constant 
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the present solution corresponds to the assumption of "constant 
properties" ) . 

For the special case of Pr = 1, equations (26) reduce merely to 

(gi)p r _i = 1 - (^2 ) px*=i “ 

Equation (25) reduces to 

^ h * ) Pr=i = 1 + fe ' 1 ) (1 ‘ + Ife "* (l 

_ f h v q U e 2 \ u 2 Ue2 

" he *Vh? “ ~ 2hJ ’ * 2h e 

which is the well-known Crocco integral for viscous flow of a gas w ith 
Prandtl number of unity. 



Velocity, Enthalpy, and Temperature Profiles 
in Physical Coordinates 


As noted previously, the relation between £ and y is determined 
by the equation 



(13) 


This equation involves T*, whereas the solution to the energy differen 
tial equation involves h*. The h(T) relationship for any gas obeying 
the equation of state p = pRT is 


h = c P a T + ^int 


(28) 


where hq n q 


■/ 


'Pint 


dT is the enthalpy for internal motion within the 


molecules, cp = 7 a R/(7 a - l) is the constant specific heat for the active 
degrees of freedom (translation and rotation), and Cp is the 


temperature -dependent specific-heat contribution for the internal degrees 
of freedom (vibration and electronic excitation). A diatomic gas, or a 
gas with linear molecules, corresponds to 7a = 7/5; & polyatomic gas 
corresponds to 7 a - 4/3- The total specific heat cp is obtained by 
differentiating equation (28) . 
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7a- 1 


R + 


R 4- ( c-n 4 Ct> 

7 a - 1 \ *vib -^elec 


For a given gas, h irrt i s a known function of T, which, in turn, is a 
known function of h. The hint( T ) relationship is different for each 
gas having internal energy within a molecule. In terms of the function 
(hint) = hint/he "the dimensionless temperature can be expressed as 


T * = -±- = 
* T~ 


°Pa Te 


hx, - hj 


which can be combined with equations ( 13 ) and ( 15 ) to yield 



where the Prandtl number independent function i}(£) is the same as pre- 
viously defined by equation ( 15 ) (tabulated in table i) , and the supple' 
mentary Prandtl number dependent functions and rj 2 (£) are defined 
by the equations 



Values of and r\ 2 are tabulated in table III. Equation (31) enables 
the velocity and enthalpy profiles in £ coordinates to be converted 
into the (x,y) physical coordinates. The enthalpy profiles together with 
the h(T) function for a given gas enable the temperature profile to be 
determined. The velocity and temperature profiles, of course, determine 
the Mach number profiles. 
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For the special case of gases having a constant specific heat (which 
would he the case of gases having no vibration within molecules , or of 
gases over the special temperature range where the vibrational energy is 
fully excited) the equations relating y and £ simplify considerably. 

The integral term in equation ( 31 ) vanishes since hj_ n t is 0. Moreover, 
h = [7/(7-1)3RT for this special case. Consequently 




and the Mach number profile is determined by 


M = Me 


u. 


Me u * 


The total enthalpy profile 


, t 1 Ty 'N 7 - 1 „ 2 

1 + l .f£ * 1 J Sl + ~ 2 ~ ^ ® 2 


(33) 


(3*0 


h _ 2 

h t + 2h e U * 

" UgJ. 


2h, 


becomes for this special case 


1 + 


^ 

Tte hte 


(I- 1 ) 


7 " 1 2 / 2 \ 

Si + — — “e ( u * + 82) 


7-1 p 

1 + — -Me a 


(35) 


Heat Transfer and Recovery Factor at Wall 


Inasmuch as the solution to the differential equations of viscous 
flow pertain to the motion only within the thin mixing layer, which is 
located a considerable distance from the wall surface, it is necessary 
to relate the rate of heat transfer at the wall to the properties of the 
mixing layer. This is done by means- of the law of energy conservation. 
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Within any closed control surface A, the conservation of energy requires 
that 


heat 
input 
across A 


increase in 
internal 
energy 
and 

kinetic 

energy 


+ work done 
by gas 
through 
normal 
pressures 
on A 


+ work done 
by gas 
through 
tangential 
stresses 
. on A 


The work contribution of tangential stresses often is not considered 
in integral formulations of the energy conservation law, but it must be 
included here because of the type of contour selected subsequently. For 
steady flow, the energy law can be converted to the form 


“ ■ //( e + f) pVndA + II ^ + II TVtdA (36) 

A A A A 

* 

where the various symbols are: 

n coordinate normal to the surface A (positive if directed outward) 

A coefficient of heat conduction (cpp/Pr) 
e internal energy per unit mass (h - p/p) 

V n velocity component normal to A (positive if directed outward) 

V-j. velocity component tangential to A (positive if shear 
stress represents work done by gas within A) 

t shear stress tangential to A (|p chji/c)y[ ) 


Along a streamline V n is 0, and along a wall V-^ is 0; hence by 
applying this conservation theorem to the closed contour formed by the 
dividing streamline SR (designated by subscript o, as sketched in 
fig. l) together with the wall boundary RS, there results 



dx-Q w = 0 + 0 + 


/7i; 


(-uo)dx 


(36a) 


where b is the width of the two-dimensional flow, and where Q w is 
the rate of heat addition to the wall (-Q w is the rate -of heat addition 
to the gas). The tangential velocity u Q has a negative sign prefixed 
inasmuch as the shear stress t does work on the gas. By solving 
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for Q w , substituting A - c p |j./Pr, and employing £ as an independent 
variable, there results 


Qw _ v r l * o /<ih* &\ 

1 e J Pr \dC 8y7 

O o 


dx 



( du* d£\ 

A d£ 


u Q dx 


(37) 


According to equation (22), G(o) = 0, hence differentiation of equa- 
tion (21) yields 


an.3 „ ^ u *o u *o"> 

= " 1 aT - 


Pr 


(38) 


Also, from equation (12) 


= p ° u ° 
‘^ y 4 v e u e xC 


(39) 


so that the integrals in equation (37) become (^since P*M-* = C = P* 0 ^* 0 ^) 


% 

b 


h.0 J PgU e P e 7-C 


-C- 


(^U-frpUttQ ' ) 

Pr 


Pr 


+ 2 


he 


U ' 

*0 *0 


(4o) 


This is a relatively simple equation inasmuch as u* o (=u Q /u e ) and 
u* o ’ ( = du*/d£ at £ = 0) are constants independent of Pr, having the 
values O.587 and 0-34l, respectively. 

Under adiabatic conditions at the wall Q, w = 0, and the value of the 
enthalpy at the wall is h aw , so that by substituting C x from equa- 
tion (24) into equation (4o) 


. X1 aw 


= 1 + 


u. 


2hc 


F 2^d) + 


4Pru 


*0 u *o 


(Sox ’ 1 


Pr 


from which it is clear that the recovery factor r 

Isci'- . iszh. mr . . iai 

-h e ( 4u * 0 u *o’ 


Fi(Cd) 


is as follows : 


t 



(M) 


(k2) 
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The dimensionless, average, heat-transfer rate per unit area q w = Q w /b7 
becomes 


p e u e (h a v “ ^-w) 


r-=- (4U* u* ' ) ?r 
/ ^ *0 *0 ' 

J Re Prp i (^ d ) 


(*3) 


Numerical computations of r from equation (42) are tabulated in 
column (a) of table IV. Values of q w i£IIZ£ computed from 

PeMhav-^) 

equation (43) are tabulated in column (c) of table IV. 


Heat Transfer When Gas Is Injected 
Into Separated Region 


It is surprisingly easy to generalize the preceding analysis in 
order to consider conditions where gas is injected into the dead-air 
region. When a mass flow m^ is injected, there no longer is an equality 
between the mass flux drawn out along the mixing layer and the mass flux 
reversed near the reattachment zone. Instead, in a steady flow, the same 
amount of mass flow m^ as is in'jected into the dead-air region also 
would pass downstream near reattachment between the dividing streamline 
and some streamline displaced at a distance y-j_ below the dividing 
streamline . 



As indicated in the sketch, where only those streamlines below the dividing 
streamline are drawn, the mass -flow parameter corresponding to the par- 
ticular streamline displaced a distance y^ from the dividing streamline 
is designated as £-j_. The equation for the balance between mass flow 
injected (near B) and mass flow escaping (near R) is • 


/ o no 

pu dy= »/p e u e |i e ZC / a? = -(^PeUgHglC 


(kk) 
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The law of conservation of energy for the contour RBSAR in the sketch 
is (see eq. ( 36 )) 



which becomes, after dimensionless variables are introduced and the 
assumption is made that the enthalpy of the injected gas (hj_) is a con- 
stant equal to the enthalpy at the wall , 


Qw 

il e' s /pe u e^e ^ ^ 




(*5) 


This equation perhaps could have been written down directly by considering 
the conservation of energy for a contour comprising the streamline BR 
and the wall RB. It also can be expressed in terms of q w = Q^/bZ as 




■t (■ 


h e «s/p e u e M-e ZC £ 


u, 


h * - h w* + 2h U* 


e - 2 >a£ 
A / 
e 


or, by using equations (21), ( 23 ), and (24) 


(46) 


q Z 


he^e^e^e^ 


= -C 1 F 1 (^ d )l 1 (^ i ) + 7 2 (^i) 


(^7) 


where the functions 1 ± and I 2 are defined as 

, £ ' r* t, / \ n A ^ 


i l(t) ,f [i . 


= f [i - Sl G)]a£ 
Fx(y j \ 


(48) 


1 2 (t) = f Cf 2 (0 - F 2 (£ a ) + u* 2 ]dt = -Ix(5)F 2 (5 a ) + f [g 2 «) + u* 2 ]a£ 

(49) 
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By substituting into equation (46) the expression for C x from equa- 
tion (24) and dividing by the rate of heat transfer without injection, 
there results 



The value for the recovery factor is 


r = F 2 (Ca) 


I 2 (£i) 
+ Ii(£i) 


(51) 


These equations depend only on the mass flow injected into the dead-air 
region, as determined by the dimensionless variable The calculation 

of the recovery factor and heat transfer in a given case requires a 
knowledge of the functions 1 ± and I 2 which are tabulated in table V, 
and the constant F 2 (£a) which depends only on the Prandtl number as 
tabulated earlier. 


Skin Friction 


By applying the momentum conservation law to the same closed contour 
as that to which the energy conservation law was applied (contour composed 
of dividing streamline SB in fig. 1 together with wall boundary BS) , 
the effective skin friction can be estimated. Inasmuch as there is no 
mass flux in or out of this contour, momentum conservation requires that 


force on wall = p x = b f R \(p. 

in x direction J \ ay J 

o o 


+ 


(p e - p)cos(n,x) 


dx 


(52) 


where (n,x) is the angle between the x direction and the outward normal 
direction, and x^ is the distance to the reattachment point. It is 
seen that the term involving pressure does not vanish as it did in the 
energy conservation integral. This is because p in the momentum inte- 
gral is multiplied by the nonzero term cos(n,x), whereas p in the 
energy integral is multiplied by the normal velocity component V n which 
vanishes for the particular contour selected. The term involving cos(n,x) 
depends on the exact shape of streamlines near reattachment , and would be 
different for each flow. Consequently, it does not seem possible to make 
any more than a rough estimate of the skin friction for the general case. 
Such an estimate can be made by disregarding the term (p e - p)cos(n,x) 
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which is zero everywhere except near the reattachment zone (where locally 
p departs from p e ) . Thus, 


C F s 


|p^ie 2 (bl) 



(53) 


By using equation (39) and noting that (p„u ) = C, the integral becomes 

* * o 





( 54 ) 


It will be seen shortly that this estimate of skin friction is the same 
as that which would be obtained by arbitrarily applying Reynolds analogy 
to the heat-transfer rate calculated for Pr = 1. 


If gas is injected into the dead-air region either at low velocity 
or in the direction normal to the mixing layer, then only one additional 
term 



must be added to equation (52) for skin friction. Conservation of momen- 
tum within the portion of mixing layer below the dividing streamline, 
however, requires that 




o 

pu 2 dy 


-00 


From equation (12) 


it is seen that 


f pu 2 dy = u e Pe"/v e Ue zc J u * d £ 


(55) 


so that the estimate of skin friction with injection can be written as 
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For the special case of no injection, integration of the u*(£) function 
(table I) yields 


(C F )^ =0 ~ 0.796 

in agreement with equation (5*0 • 

RESULTS AND DISCUSSION 

Comparison With Results for Laminar Boundary Layers 


Recovery factor.- A good way to visualize the results of the pre- 
ceding analysis is to compare them with corresponding calculations for 
a flat plate over which a constant-pressure laminar boundary layer flows . 
The two functions g-^ and g 2 , which combine linearly to yield the 
enthalpy-velocity relationship for a separated laminar mixing layer 
(according to eq. ( 25 )) are shown in figure 2. Also shown for purposes 
of comparison are the corresponding two functions for a laminar boundary 
layer (according to the analysis of refs. 5 and 8). For Pr = 1 the 
g 1 and g 2 functions are equal to 1 - % and u*(l-u*), respectively; 
these quantities also are shown in figure 2. As might have been antici- 
pated for Pr near 1, the component functions for enthalpy are much the 
same for the two types of viscous layer. Since the recovery factor is 
merely a measure of the enthalpy when the velocity is zero, it too should 
be nearly the same for the two types of laminar flow. The recovery factor 
for a flat plate in an incompressible stream is very nearly equal to n/Pt, 
as first calculated by Pohlhausen (ref. 9 ). In a compressible stream the 
recovery factor is unaffected by variation of Mach number within the 
framework of assumptions made in the present analysis. A comparison of 
4¥r with the values of recovery factor computed herein for separated 
laminar flows is made in the following table: 


Recovery factor in 
separated flow 

*/pr 

Pr 

r 

0.1 

0.361 

0.316 

.25 

• 504 

.500 

.5 

.712 

• 707 

.72 

.850 

.849 

1.0 

1.000 

1.000 

1.5 

1.228 - 

1.225 

2.0 

1.423 

1.414 

10.0 

3.27 

3.16 
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It is seen that the recovery factor for pure laminar separated regions 
is approximately equal to VPr. Consequently the recovery factors for a 
separated mixing layer and an attached boundary layer are essentially the 
same . 

Rate of heat transfer .- A comparison of the corresponding rates of 
heat transfer yields different results. For the flow of a compressible 
laminar boundary layer over a flat plate under the same framework of 
assumptions as employed herein, the average rate of heat transfer (q^)-^ 
according to references 5 and 8 is < 

/*E = Yq, ^ Q ) i (x 0.664 Pr‘ 2/3 for Pr near 1^) (57) 

Pe u e( h w- h awW c Pr V ' 


A comparison with the rate of heat transfer for the pure laminar type of 
separated flow is made in the following table; 


Separated mixing 
layer 

Boundary layer 

3w 


q v / 0 

(q„) b ^Re/c 

Pr 


p e u e (h w - h aw ) 

P e u e (^w “ ^aw) 

0.1 

0.833 

2.70 a 

0.31 

.25 

.674 

1.57 a 

•43 

-5 

.526 

1.03 

*51 

.72 

.457 

.820 

• 56 

1.0 

*399 

.664 

.60 

1.5 

*335 

.510 

.66 

2.0 

.293 

.421 

*70 

10.0 

.135 

.146 

.92 


a The calculations of boundary -layer flow for 


small Prandtl numbers (e.g., for Pr = 0.1 
and 0.25) are based primarily on appro- 
priate small Pr expansion formulae developed 
by P. Lagerstrom and H. Liepmann in some 
unpublished research. 

In calculating the ratio q w / (q^-)-^ ^ is assumed that the values 

of p e , u e , Re, and h w /h~ are the same for the two cases. It is to be 
noted that the ratio qw/(q w ) b ^ of the two heat fluxes is independent 
of Mg, Re, and h w /h e . For air (Pr « 0.72), the calculated heat transfer 
in a separated laminar mixing layer is 0-5 6 of that for a corresponding 
attached laminar boundary layer. 
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Attention is called to the fact that the use of the integral form 
of the energy conservation law enabled only the average (or over-all) 
rate of heat transfer q^. to be calculated and not the local rate. In 
this respect, the present analysis of separated flows is unlike the 
analysis of an attached laminar boundary layer wherein both the local 
distribution and the average rate of heat transfer can be calculated. 

The local distribution of heat flux along the wall of a separated flow 
would depend on the details of the vortex-like motion within the region 
of dead air lying between the wall and the thin mixing layer. Without 
considering these details together with the particular shape of the wall, 
the theory cannot provide information about the local rate of heat 
transfer . 

Skin friction .- Since (Cp)-^ for an attached laminar boundary layer 

is 1 . 328 C /Re (see ref. 5), it follows from equation (5*0' that the esti- 
mated value for effective friction of a separated laminar mixing layer is 
0.80/1.328 = 0.60 of that for the corresponding attached laminar boundary 
layer. At Pr= 1 the heat flux, as tabulated above, also is 0.60 of that 
for the corresponding attached laminar boundary. Thus the estimate of 
skin friction is the same as that which would be obtained by arbitrarily 
applying Reynolds analogy. 

Velocity and temperature profiles .- To illustrate further the rela- 
tive characteristics of separated layers and attached boundary layers , a 
comparison of both types of layer for Me = 10 is presented in figures 3 
and 4 showing velocity and temperature profiles, respectively. The two 
extreme conditions of Tv/T e = 1 and Tv = Taw are considered, as well as 
a third condition representing T v /T e = 4. These examples illustrate both 
some differences and similarities between a separated mixing layer and an 
attached boundary layer. For example, mixing layers are several times 
thicker t han boundary layers, yet the maximum temperature attained within 
each layer for cold-wall conditions is essentially the same. 

Flow with gas injection .- Numerical results illustrating the effect 
of mass injection on recovery factor in separated flow are presented in 
figure 5 for the values Pr = 0.5* 0.72, and 1.0. It is evident that 
for Pr = 1 the recovery factor is unaffected by mass injection and is the 
same, for the separated mixing layer as it is for the attached laminar 
layer. When Pr < 1, however, the recovery factor is lowered substantially 
bty- injection, for both mixing layer and boundary layer. The curve repre- 
senting the laminar boundary layer with mass injection (dotted line in 
fig. 5) is taken from the calculations of Low (ref . 10) which are based on 
the same viscosity- temperature relationship as is used herein. 

The effect of gas injection on heat-transfer rate can be conveniently 
illustrated through consideration of several special cases which simplify 
the general equation (50) . First, for the special case of h v /h e = 1 
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(which would imply T v /T e = 1 if the specific heat were constant), then 
equation (50) becomes 


Qw 

(Qw)^ i= o 


bw~b e 


F a (C d )I 1 (0) + I 2 (0) 


(58a) 


which is independent of Mach number. Also, if the Mach number is very 
large so that Mg 2 > > (hw/he - l) (or ■ Ue 2 /2he > > (hw/he - l))> then 
equation (50) again reduces to the same expression 


Qw 


(Qw) 


£i=o 


Me=< 


Fa(Cd)li(Si) + la(Cl) 
F 2 (C d )I 1 (0) + I 2 (0) 


(58b) 


For gases, expressions (58) also are nearly independent of the Prandtl 
number, as the curves in figure 6 illustrate for Pr = 0.5> 0.72, and 1.0. 
-If a different limiting case is considered, namely, Me~> 0, then 
u e 2 /2h e < < (b w /h e - 1), and equation (50) becomes 




Me=o 


ii(Si) 

ii(o) 


(59) 


Numerical calculations show this ratio to be only slightly higher than 
the ratio given by equation (5 8)* For example, at -£i = 0.4, the ratio 
in equation (58) representing Mg- — >°°, or h v = he , is 0.49 (for Pr = 0.72), 
whereas the ratio in equation (59) representing Me — >0 is 0.52. For 
practical purposes, then, the curves in figure 6 are applicable to a 
wide variety of conditions . 

As perhaps might be anticipated, the effect of mass injection on the 
esti ma ted skin friction is much the same as the effect on heat transfer. 

At Pr = 1, for example, g x is equal to 1 - u*, as previously noted, so 
the combination of equations (48) and (56) yields for the ratio of estimated 
skin friction with injection to that without injection 


Cp 


( c f)^ 1=( 


Pr=i 


ii(Si) 

li(0) 


(60) 


This ratio is identical to the corresponding ratio of heat-transfer rates, 
as indicated by equation (59) for low-speed flow. Inasmuch as the ratio 
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in equation (59) is practically the same as that in equations (58) * it 
follows that the effect of mass injection on estimated skin friction is - 
for practical purposes - essentially the same as the effect on heat- 
transfer rate (as represented by the curves in figure 6) . 


Computational Checks and Extensions of Analysis 


It was found possible to obtain an independent check on the internal 
consistency of the numerical calculations by employing certain integral 
considerations other than those used in the main analysis. An example of 
such a check already has been observed. The quantity CjWRe/C for sepa- 
rated flow without injection is equal to. the differential expression 
u* 1 , according to equation (5*0 > ye t also should be equal to the 

o 

u*d£, according to equation (56). The table 

of u*^) yields values of 0.801 and 0.796, respectively, for these two 
independent expressions. This is considered to be adequate agreement for 
present purposes. Similarly, independent integral checks (usually to 
within a few-tenths of 1 percent) of the differential expressions for both 
recovery factor and heat-transfer rate were obtained for all Prandtl num- 
bers. Details are presented in Appendix A. 

Probably the most significant feature of the integral method employed 
to check computations is that it yields equations which also can be applied 
to turbulent separated flows, provided the velocity profiles in the turbu- 
lent mixing layer are known. Details of the application to turbulent 
separation are presented in Appendix B. It will suffice here to note that 
velocity-profile data for turbulent mixing layers at high Mach numbers are 
not yet available, so the numerical calculations are restricted to low 
Mach numbers where such data are available. The calculated results for 
low- speed flow indicate that the heat -transfer rate to a separated turbu- 
lent flow is much higher than that to a comparable attached turbulent 
boundary layer. This result for turbulent flows contrasts sharply with 
the corresponding result for laminar flows. 

Two other extensions can be made to the analysis. One is to axially 
symmetric flow, which is a rather simple extension and is presented in 
Appendix C. The other is to separated flows wherein the boundary-layer 
thickness at separation &s is sizable. This is not a simple extension. 
Equation (36a) would apply directly to this more general case, but great 
difficulty would be encountered in computing (dT/8y) Q , (8u/8y) 0 , and u Q 
at various stations along the dividing streamline of the mixing layer. 

In effect, the absence of similar profiles for the case where 5 S . is not 
zero would require that partial differential equations be solved for 
u(x,y) and T(x,y), rather than ordinary differential equations. 


integral expression 2 
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Some Practical Aspects of Gas Injection 


The practical aspects of gas injection for separated flows are 
expected to be quite different than for flows of an attached laminar 
boundary layer. It is well known that the injection of gas through a 
porous wall unfortunately has a strong destabilizing effect on the lami- 
nar boundary-layer flow along that wall. There is no evident reason, 
however, to expect that the gas injection into a separated region through 
a porous wall, which is displaced a considerable distance from the sepa- 
rated laminar mixing layer, would necessarily have a strong destabilizing 
effect on the stability of the mixing layer. Such injection, if properly 
done, might have a favorable effect on stability inasmuch as the pressure 
rise near reattachment is reduced when gas is injected, and this pressure 
rise is important to the stability of the laminar flow. Some experiments 
on the effect of injection on stability clearly are in order. 

If a mass flux of gas equal in magnitude to 


pu dy = J p e u e pe ZC J d£ = J p e u e p e ZC (65) 

£d 

is injected into a separated region, then no air would be reversed near 
the reattachment zone . For this particular quantity of injection the 
heat transfer (and the estimated skin friction) would be zero. In mathe- 
matical terms, for = £ d , then Ii(£ d ) = la(£ d ) = 0 an(i Qw there- 
fore zero (see eq. (50)). It is interesting to consider an example in 
order to obtain a physical feeling for the magnitude which this particular 
rate of mass injection would represent in a practical ease. Consider 
first, as sketched, a cone with 30° total included angle, a base radius 




of 1 foot, flying at a Mach number of 20 at an altitude of 100,000 feet, 
and containing initially l/3 of its volume as stored liquid air (or liquid 
nitrogen). For the case = £ d , for which Qw ~ 0, the cone would 
travel in steady flight a distance of approximately 1.4x10® body diameters 
(or about 500 miles) before the stored liquid had been completely injected 
into the separated region. Along this length of flight path, the heat 
transfer would be zero. After all mass had been expended, the rate of 
heat transfer would be O.56 of that for a cone of 15° semiangle (provided 
the separated layer remained laminar). 
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Inasmuch as the quantities of mass injection required are relatively 
low, the scheme of mass injection into separated regions may have practi- 
cal application. The degree of practicality depends upon the Reynolds 
numbers up to which separated mixing layers can be maintained laminar . 
Experiments are required to ascertain this under the conditions of hyper- 
sonic flight where the wall is cool compared to stagnation conditions and 
where gas is injected into the separated region. 


CONCLUSIONS 


The conclusions which follow have been obtained from an analysis of 
separated laminar flows wherein the thickness of the boundary layer is 
zero at the separation point: 

1. The recovery factor in regions of separated laminar flow is 
approximately equal to the square root of the Prandtl number (to within 
1 percent over the range of Prandtl numbers between 0.25 and 2.0) and, 
hence, is essentially the same as that of an attached laminar boundary- 
layer flow. 

2. The calculated rate of heat transfer from a separated laminar 
mixing layer is less than that from an attached laminar boundary layer at 
corresponding values of Mach number, Reynolds number, and wall to stream 
temperature ratio; the ratio of the heat flux in the separated flow to 
that in the attached flow is a function only of the Prandtl number, and 
has the numerical value of O .56 for Pr = 0»72 (the value for air). 

3 • Injection of gas into the dead-air region of a separated laminar 
flow is calculated to have a powerful effect in reducing heat -transfer 
rates, inasmuch as a moderate quantity of injection can reduce the heat 
flux to zero. 

A. The effective skin friction for separated laminar flows is 
estimated to be about 0.6 of that for an attached laminar boundary-layer 
flow, and is affected by gas injection in essentially the same way as 
heat transfer is . 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
June 25, 1956 
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APPENDIX A 

CONTOUR INTEGRAL CHECKS ON ACCURACY 
OF NUMERICAL CALCULATIONS 


A check for internal consistency on the numerical calculation of 
u*(£) can he obtained by noting that the conservation lav of linear 
momentum in the x direction , when applied to a contour enclosing the 
mixing layer, requires that 


J pu 2 dy = J pu(u e - u)dy 


In (x,£) coordinate system, this requirement becomes 


f u*d£ = /^ (l - u*)d£ 
1 


(Al) 


Numerical evaluation of these two integrals from the u*(£) solution 
yields ' 


J u*d£ = 0.398 

£d 


(A2) 


and 



(A3) 


These two integrals agree within 0.29 percent. As a further check on 
u^, it may be noted from the differential equation (ll) that when u* = 0, 
(du*/d£)^ must be the same as -£3/2. The values in table I show that 

(du*/d£)£ = +0.620 and -£3/2 = +O.617. These values check within 
0.5 percent. 


Checks on the calculations of r and Q w can be obtained from the 
energy conservation law by considering a contour enclosing only the dead 
air region. The flux of energy fed into the dead-air region near the 
reattachment zone is 


b 





pu dy 



MCA TN 3792 


29 


This energy flux is transported by the mass flux 



pu dy, which 


-00 

must be equal to the mass flux reversed by the pressure rise near the 
reattachment zone. (This mass flux, in turn, must also be equal to the 
mass flux re-entering the mixing layer after some heat energy has been 
imparted to the wall.) When this mass flux leaves the dead-air region 
and enters the mixing layer, it transports an energy flux of 


b 



l 

h wPw^w^ 


which, from the continuity equation, is equal to 


b f h^pu dy 


Consequently, the rate of heat transfer to the wall must be the difference 
between these two energy fluxes . 


Qw 


= b J (h - h v + ~“^pu dy 


In dimensionless variables. 


Qw — */pe^e^e ^ J' “ ^w* + 2hg^ ^ 

^d 


(A4) 


(A5) 


By substituting h* from equation (21) , and setting = there 

is obtained 


*w 


/Re 

p e u e h e 7 c 



[F 1 (U-F i a d )]a$ + ^- J 

a 


[f 2 (0 -F a (C d )]at + 


(a6) 



30 


NACA TN 3792 


Under adiabatic -wall conditions q v = 0, and h w - h aw ; hence by substi- 
tuting C ± from equation (24) 


^aw _ -j_ + Ue 
hp 2h e 


3d£ 


it p 2 (5) - F 2 (5 d ) +'u» 

F 2 U d ) - F i(?d)— 1 T 75 / 

J [Fi(5) - Fi(£d>^ J 


from which it follows that an expression alternate to equation (42) for 
the recovery factor is 


s: 


, [f 2 (0 - F a (5 d ) + u a ]d£ (Q) 

r - F a (£ d ) - Fl (C d ) -Ji™ * = F a (6 d ) + (AT) 


/ [F x (5) - F 1 (5 d )dE 


ii(o) 


where the functions I 2 and I x are, as defined previously by equa- 
tions (48) and (49) 

.i 


: X (E) = f [i - 6l (£)K 


i 2 (0 = r [f 2 (£) - F 2 (£ d ) + u* 2 ]d£ 

ta 

Expression (A7) should be numerically equal to the different expression 
of equation (42) . The corresponding alternate expression for the rate 
of heat transfer is 


<lw 


P e u e(haw - hw) 


- 71 7 11 ■ “ <!)W 


(A8) 


which should be equal to the expression of equation (43) . 

Values of r computed from the alternate equation (A7) are tabulated 
in column (b) of table IV. Compared with the values in column (a) origi- 
nally computed from equation (42), it is seen that the alternate inde- 
pendent calculations agree rather well with the original computations . 

The largest discrepancy amounts to 0.2 percent for values of Prandtl num- 
ber less than 2. For a Prandtl number of 10, the two computations differ 
by 3 percent. In this latter case, the values from the integral method 
(eq. (A7) ) are regarded as less accurate, inasmuch as their evaluation 
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Pr* 

involved, several quadratures of functions such as u* which, for 
Pr = 10, vary so rapidly that high accuracy could not be obtained with 
the intervals selected. Nevertheless, the over-all agreement between the 
two independent methods of computation is regarded as satisfactory. 
Similarly, a satisfactory check on numerical computations of q w by the 
two independent methods is obtained, as evidenced by comparison of col- 
umns (c) and (d) in table IV. Compared with the original calculations 
from equation (lj-3) (column (c)), the calculations from the alternate 
equation (A 8) (column (d)), agree to within 0*5 percent for all values of 
Prandtl number. 
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APPENDIX B 

HEAT TRANSFER IN SEPARATED TURBULENT FLOWS 


Although complete calculations cannot at present he made for sepa- 
rated, compressible, turbulent mixing layers, certain ol the ideas 
developed nevertheless can be applied to such cases. For this purpose 
it is assumed that the Prandtl number is unity. Equation (A^) was 
developed from general energy considerations and is applicable also to 
turbulent flows . 



-oo 


hw 


+ 



pu dy 


(Bl) 


The assumption Pr = 1 enables the relation between enthalpy and velocity 
to be expressed as the Crocco integral, 


** = 1 + Gr - Mt 1 ■ + u *( x - u *) 


= 5a + uii + £ 

h P *\ 2h e 


- h w ) - u 


2 

' 2hf 


(B2) 


so that the average rate of heat transfer becomes, after some algebraic 
manipulation wherein the relationship hav* = (u e 2 /2h e ) + 1 is used and 
the dimensionless variable | s dy/x is introduced, 


*w 


Pe u e(^aw " ^w) ^ 


p* u * Sd6 


1 


u^di 


ah 


w. 


* - 


(B3) 


- U, 


2hw 


where the constant a is inversely proportional to the rate of spread 
of the turbulent mixing layer. Clearly, it is necessary to know the 
velocity distribution u*(£) within a turbulent mixing layer, and the 
rate of spread a, before the average rate of heat transfer can be cal- 
culated. Such knowledge, unfortunately, is not yet available except for 
two cases . 

For the limiting case of zero Mach number, the velocity distribution 
and rate of spread are known from the paper of Tol'lmien (ref. 11). By 
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substituting Tollmien’s velocity distribution (cr = 12 at Re = 2.5x10 s 
from his experiments) into equation (B3) there results 


— at = 0.012 

Pe u e (h a w - ^w) 

which leads to an interesting result when comparison is made to a corre- 
sponding attached turbulent boundary layer . At the same Reynolds number 
as in Tollmien’s experiments, the corresponding heat-transfer parameter 
for a turbulent boundary layer of constant pressure would be approximately 


‘ 5w 

-Pe u e(* 1 a'w 


h v ) 




0.0019 


which indicates the heat transfer in separated turbulent flow not to be 
smaller than in a turbulent boundary layer, but to be, in fact, 

— = 6.3 times as large. This comparison is in sharp contradistinc- 
0.0019 

tion to the analogous comparison for the laminar case where the corre- 
sponding ratio is about 0.6 rather than 6 . 3 . 


It is to be noted that the data available at present for compressible 
turbulent mixing layers indicate the Mach number to have a pronounced 
effect on the integral in equation (B3) . The data of Gooderum, Wood, 
and Brevoort (ref. 12) at a Mach number of 1.6 indicate that the rate of 
spread (about 9° angle) of turbulent mixing layer is much less than at 
low speed (about l4° angle); hence the integral in equation (B3) also 
would be proportionately less. Consideration also of the density change 
on going from M = 0 to M = 1.6 would then yield a value for this integral 
57 percent less than that at low speed, which is equivalent to a heat 
flux in turbulent separated flow of about 2.8 times that in a correspond- 
ing turbulent boundary layer. The ratio of heat fluxes, amounting to 
6.3 at Mg = 0 , and 2.8 at Mg = 1 . 6 , obviously is strongly dependent on 
Mach number. Consequently, if the marked trend persists, a separated 
turbulent flow might not have a greater heat flux than a turbulent boundary- 
layer flow at sufficiently high Mach numbers . Experiments on the rate of 
spread and on the velocity distribution within turbulent mixing layers at 
high supersonic speeds are required before this can be ascertained. 
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APPENDIX C 


APPLICATION TO AXIALLY SYMMETRIC FLOW 


Although the analysis as developed applies directly to two- 
dimensional flow only, it can he easily extended to axially symmetric 
flow by employing Mang ier 1 s transformation. This transformation for a 
c onstant -pres sure , separated, laminar mixing layer is identical to that 
for a laminar boundary layer. For a given Reynolds number and Mach nuim 
ber, the mix ing-layer thickness in axially symmetric flow is I/nTT times 
that in two-dimensional flow. Consequently, the heat-transfer rate, in 
place of equation (V?) > is given by the equation 


Qw 

he •'J Pe u eM , e^G 


2 * r b 

4T 


Ei 


/(■ 


h* ~ hw* + 


u e 2 

2he 



(Cl) 


and hence the average heat -transfer rate per unit area 


q.w s QwAV is 



he Vp e ^e^e zc 




(C2) 


which is seen, by comparison with equation (^-6) to be 2-/*J~3 times as large 
as for two-dimensional flow, just as it is in the case of comparison of 
an attached boundary layer on a cone with the corresponding laminar bound- 
ary layer on a plate. Consequently, the ratio q w /(q w )^^ tabulated pre- 
viously for two-dimensional flow also applies directly to axially symmetric 
flow. If gas is injected into an axially symmetric region of separated 
flow, then for the same fractional reduction in heat-transfer rate (same 

(2//3)(m i /jtr b z) 

value of q.w/(q. w )£ i=0 ) the value of = — u n C/l ^ aXially 

m ± /hl 

symmetric flow would be the same as the value of . in 

Vp e u e p e C/l 

two-dimensional flow. 
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TABLE I. - VELOCITY FUNCTION u *(£) AND y DISTANCE FUNCTION T)(£) 
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TABLE II.- ENTHALPY-DISTRIBUTION FUNCTIONS 
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TABLE II.- ENTHALPY-DISTRIBUTION FUNCTIONS - Concluded 
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TABLE III.- SUPPLEMENTARY y DISTANCE FUNCTIONS - Concluded 
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TABLE IV.- COMPARISON OF INDEPENDENT METHODS OF CALCULATING 
RECOVERY FACTOR AND HEAT-TRANSFER RATE 


Pr 

Recovery Factor 

Heat Transfer 

(cj 

q^Re/c 

(d) 

qw^fee/C 

(a) 

r from eq. (42) 
4Pru^. 2 u* 1 

-MSd)- 2 — 

( ltu *o u *o') 

(*>) 

r from eq. (AT) 
I 2 (0) 

= ^ ( ^ ) + xI(0) 

Pe u e(baw “ 
from eq. (43) 

K u *o'f 

PrFi(?d) 

Pe u e(baw 
from eq. (A8) 

= i x (o) 

0.1 

0.361 

0.360 

0.833 

O .836 

.25- 

.504 

.503 

.674 

.674 

• 5 

.712 

.713 

.527 

.525 

.72 

.849 

.849 

.458 

.456 

1.0 

1.000 

1.000 

•399 

.398 

1.5 

1.228 

1.227 

.334 

.334 

2.0 

1.424 

1.422 

.293 

.293 

10.0 

3.27 

3.40 

.135 

.136 


TABLE V.- FUNCTIONS APPEARING IN EQUATIONS FOR GAS INJECTION 


h 



Pr = 

0-5 

Pr = 

0.72 

Pr = 

1.0 



ii 

I 2 



- 1.233 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

- 1.2 

.002 

.002 

.001 

.001 

.000 

.000 

- 1.1 

.019 

.013 

.011 

.010 

.006 

.006 

- 1.0 

.044 

.031 

.028 

.025 

.017 

.019 

-.9 

.076 

•055 

.051 

.047 

.033 

.039 

-.8 

.112 

.083 

.080 

.074 

.055 

.065 

-.7 

.152 

.116 

.113 

.106 

.082 

.097 

- .6 

.196 

•153 

.151 

.144 

.115 

.134 

-.5 

.244 

.194 

.193 

.186 

.151 

.178 

-.4 

.295 

.241 

.239 

•233 

•193 

.226 

-•3 

.349 

.291 

.288 

.285 

.238 

.279 

-.2 

.405 

.346 

.341 

• 341 

.288 

.337 

-.1 

.464 

.405 

.397 

.402 

• 341 

.400 

0 

.525 

.468 

.456 

.467 

.398 

.467 
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Figure 3, -Velocity profiles in separated laminar mixing layer and in laminar boundary 
layer for various wall temperatures; M,*IO, P r *0.72. 
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Figure 5.- Effect of mass injection on recovery factor. 



Figure 6.— Effect of mass injection on heat transfer for case h w /h e = l 
(also applies to case Me-**® for arbitrary h*). 


NACA - Langley Field, V*. 





